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ABSTRACT 



Aims. We update the constraints on the time variation of the fine structure constant a and the electron mass m e , using 
the latest CMB data, including the 7-yr release of WMAP. 

Methods. We made statistical analyses of the variation of each one of the constants and of their joint variation, together 
with the basic set of cosmological parameters. We used a modified version of CAMB and COSMOMC to account for 
these possible variations. 

Results. We present bounds on the variation of the constants for different data sets, and show how results depend on 
them. When using the latest CMB data plus the power spectrum from Sloan Digital Sky Survey LRG, we find that 
a/ao = 0.986 ± 0.007 at 1-a level, when the 6 basic cosmological parameters were fitted, and only variation in a was 
allowed. The constraints in the case of variation of both constants are a/ao = 0.986±0.009 and m e /m e o = 0.999±0.035. 
In the case of only variation in m e , the bound is m e /m e o = 0.964 ± 0.025. 
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1. Introduction 

The variation of fundamental constants over cosmological 
time scales is a prediction of theories that attempt to unify 
the four interactions in nature, like string derived field the- 
ories, rela ted brane-w orld theories and Kaluza-Klein the- 
ories (see Uzan (2003) and references therein). Many ob- 
servational and experimental efforts have been made to 
put constraints on such variations. Most of the reported 
data are consistent with null variation of fundamental con- 
stants. Although there have been recent claims for time 
variation of the fine structure constan t (a) and of the pro- 
ton to electron mass ratio (/i = ^f-) (|Murphv et al.l [20031 ; 
iReinhold et al.|[2006h. independent analyses of similar data 
give null results (ISrianand et al.l 120041: iKing et al.l 120081; 
iThompson et al.l [2009; M alec et al.l l2010t ). On "the other 
hand, a recent analysis of ammonia s pectra in the Milky 
Way suggests a spatia l variation of [i (jMolaro et al.l 120091; 
iLevshakov et al.ll2009l ). 

Unifying theories predict variation of all coupling con- 
stants, being all variations related in general to the rolling 
of a scalar field. Therefore, the relationship between varia- 
tions of coupling constants depends on the unifying model. 
In this paper we adopt a phenomenological approach and 
analyse the possible variation of a and/or m e at the time 
of the formation o f neutral hydrogen witho ut assuming any 
theoretical model. Nakashi ma et al.l (|2010fl have considered 
also the variation in the proton mass (m p ). This quantity af- 



fects mainly the baryon mass density and the baryon num- 
ber density. Their results confirm the strong degeneracy 
with the baryon density. Therefore, we will not consider 
the variation in m p in this work. 

Cosmic microwave background radiation (CMB) is one 
of the most powerful tools to study the early universe and 
in particular, to put bounds on possible variations in the 
fundamental constants between early times and the present. 
Changing a or m e at recombination affects the differential 
optical depth of the photons due to Thompson scattering, 
changing therefore Thompson scattering cross section and 
the ionization fraction. The signatures on the CMB angular 
power spectrum due to varying fundamental constants are 
similar to those produced by changes in the cosmological 
parameters, i.e. changes in the relative amplitudes of the 
Doppler peaks and a shift in their positions. Moreover, an 
increment in a or m e decreases the high-i? diffusion damp- 
ing, which is due to the finite thickness of the last-scattering 
surface, and thus, increases the power on very small scales 
([Kaplinghat et al.lll999l; lHannestadlll999l ). 

Recent analysis of CMB data (earlier than the 
WMAP seven-year release) includ ing a possible varia- 
tion in a have been p erformed by Scoccola et ahl (2008, 
20091): iMenegoni et al.l (|2009f k iNakashima et al.l (|2010l) : 
Martins et ahl (I2010D . and inc luding a possib l e variation i n 
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m e have been perform ed by IScoccola et ahl (|2008l l2009h : 
INakashima et al.l (|2010l) . 

In our previous works, we have also analyzed the de- 
pendence of the updated recombination scenario (that in- 
cludes the recombina tion of helium, and was implemented 
in Recfast following lWong et ahl ([2008)) on a and ra e , and 
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show that these dependencies are not relevant for WMAP 
data. 

In this paper we adopt a phenomenological approach 
and analyse the possible variation in a and/or m e without 
assuming any theoretical model. We use WMAP seven-year 
release, together with other recent CMB data. We also com- 
bine CMB data with other cosmological data sets: i) the 
power spectrum of the Sloan Digital Sky Survery DR7 LRG, 
ii) a recent constraint of the Hubble constant Ho with data 
from the Hubble Space Telescope. In section [3] we describe 
the method and data sets we used in the statistical anal- 
ysis. We present and discuss our results in section [3] We 
conclude in section UJ 



2. Statistical Analysis 

We performed our statistical analysis by exploring the pa- 
rameter space with Monte Carlo Markov chains gener- 
ated with the CosmoMC code (iLewis fc B ridle 2002) which 
uses the Boltzmann code CAMB ([Lewis et al.l [20001 and 
Recfast to compute the CMB power spectra. We modi- 
fied them in order to include the possible variation in a and 
m e at recombination. 

We use data from the WMAP 7-year temper- 
ature and tempe rature-polarization power spectrum 
(lLarson et al. 20101). and ot her CMB ex periments such a s 
CBI ijReadhead et al.l l2004f). ACBAR kuo et al.l j 200? 



BOOME RANG (iPiacentini et al.ll20"06t Uones et al.ll2 0061. 
BICE P (|Chiang et all 12009ft and QUAD ([Brown et al.l 
I2009T) . In order to reduce degeneracies of the cosmologi- 
cal parameters, we combine the CMB data sets with other 
cosmological data: i) the power spectrum of the Sloan 
Digital Sky Survey LRG ([Reid et al.ll2009ft and ii) the re- 
cent constraint on the Hubble constant, Hp = 7 4.2 ±3.6 
km s" 1 Mpc -1 , presented by IRiess et all (|2009D . We did 
not consider supernovae t ype I data, because these data 
(see iKowalski et al.l ([2008) for example) are obtained as- 
suming that the constants have their present values at the 
time corresponding to the observing redshift. However, if at 
high redshift a has a value different than the present one, 
the light curves of the SNs could be affected by these varia- 
tions. Since there is no agreement in the scientific commu- 
nity about the claimed variations at redshifts 0.22 < z < 3, 
and the supernovae data set includes data at redshifts 
0.001 < z < 1.4, we decided not to include this data set 
in our analysis. 

We have considered a spatially-flat cosmological model 
with adiabatic density fluctuations, and the following pa- 
rameters: 



P = ( n b h 2 , ilcDMh 2 , O, r, — , — ,n s ,A. 

a m e0 



where fl^h 2 is the baryon density and ilcDAih 2 is the dark 
matter density, both in units of the critical density; O gives 
the ratio of the co-moving sound horizon at decoupling to 
the angular diameter distance to the surface of last scat- 
tering (and is related to the Huble constant Ho); t is the 
reionization optical depth; n s the scalar spectral index; and 
A s is the amplitude of the density fluctuations. 

We have performed statistical analyses using the data 
mentioned above and considering variation of only one con- 
stant (a or m e ) and variation of both constants. We present 
our results in the next section. 



3. Results and Discussion 

Results for the variation of the constants in the case 
when only one constant is allowed to vary are shown in 
Table Q] and for the case when both are allowed to vary, 
are presented in Table The obtained values are consis- 
tent with no variation of a or m e at recombination. The 
obtained err ors are at the s ame percent level than those 
obtaine d by IScoccola et al.l ([20081 I2009H : iMenegoni et al.l 
(|2009ft ; iMartins et al.l (|2010t) using WMAP-5 year release. 
The parameter space has higher dimension when both con- 
stants are allowed to vary. Therefore, limits on a and m e 
are more stringent in the case were only one constant is 
allowed to vary. Results for the cosmological parameters 
have similar values for all of the analyses. Therefore, we 
only report the values obtained in the case where both a 
and m e were allowed to vary and the data from CMB and 
the power spectrum of the SDSS DR7 were considered (see 
Table [3]) . The mean values and errors for the cosmological 
parameters are in agreement w ithin 1-g with those obtained 
by the WMAP collaboration ([Larson et all 12010(1 with no 
variation of fundamental constants. 



Table 1. Mean values and l-cr errors for the analysis with 
variation of only a, and only m e . 



Data set 


a/ao 


m e /m e o 


all CMB 


0.987±°:°oo°9 


0.9831°;°^ 


all CMB + H 


0.9981°:°°? 


i-oi2l°:°^ 


all CMB + Sloan P{k) 


0.986 ± 0.007 


0.964 ± 0.025 



Table 2. Mean values and 1-a errors for the analysis with 
the joint variation of a and m e . 



Data set 


a/ao 


m e /m e0 


all CMB 


0.986 ± 


0.010 


i.uio_ 074 


all CMB + H 


0.986 ± 


0.010 


1.044 ± 0.029 


all CMB + Sloan P(k) 


0.986 ± 


0.009 


0.999 ± 0.035 



In Fig. Q] we show the 68% and 95% c.l. constraints 
for a/ao versus Hq, for the analysis of the variation of a 
alone. The results correspond to different data sets: all the 
CMB data alone' all the CMB data plus the Ho prior taken 
from IRiess et al.l (|2009h : and all the CMB data plus the 
power spectrum f rom Sloan Digital Sky Survery DR7 LRG 
([Reid et aT]|2009h . The large degeneracy between a/a and 
Ho from CMB data is reduced when another data set is 
added. However, since the value of Ho obtained from the 
extra data sets are different, the obtained constraint on 
a/ao depends strongly on the data chosen for the analysis. 
Nevertheless, the results are consistent within l-cr. 

In Fig. [5] we present the constraints for a/ao versus 
t and in Fig. [3] we present the constraints for a/ao versus 
Q,bh 2 . There are degeneracies among these parameters. The 
contours change because of the different mean value of a/ao 
obtained with different data sets. 
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Table 3. Mean values and lcr errors for the cosmological 
parameters using all CMB data and the SDSS DR7 power 
spectrum. Hq is in units of km s _1 Mpc -1 . 



Fig. 1. 68% and 95% c.l. constraints for a/ao versus Ho, 
for the analysis of the variation of a alone. Results from 
different data sets. 



Fig. 3. 68% and 95% c.l. constraints for a/ao versus Qbh 2 , 
for the analysis of the variation of a alone. Results from 
different data sets. 
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Fig. 2. 68% and 95% c.l. constraints for a/ao versus r, 
for the analysis of the variation of a alone. Results from 
different data sets. 
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Ho is larger than between a/ao and Ho , making impossible 
to find reliable constraints using CMB data alone. When 
another data set is added, the bounds result tighter, but 
the mean value for m e /m e o depends strongly on which data 
set was added. Results are marginally consistent at 1-a. 



Fig. 4. 68% and 95% c.l. constraints for m e /m e o versus Hq, 
for the analysis of the variation of m e alone. Results from 
different data sets. 
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In Fig. 0] we present the result for the case where only m e 
was allowed to vary. The degeneracy between m e /m e o and 
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The constraints on m e /m e o versus r are shown in Fig. [SI 
and on m e /m e o versus flbh 2 are shown in Fig. [JJ] In both 
cases, the results depend on the data set added to CMB 
data in the statistical analysis. 

In Fig. [7] we show the posterior distribution for a/ao 
and m e /(m e )o, for the case of joint variation of these quan- 
tities, marginalized over the cosmological parameters. The 
results correspond to different data sets. The difference in 
the contours is mainly due to the large degeneracy of m e 
and Hq, and the different Hq values derived from the Sloan 
power spectrum and from the Hq prior. We see that the 
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Fig. 5. 68% and 95% c.l. constraints for m e /m e o versus r, 
for the analysis of the variation of m e alone. Results from 
different data sets. 
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Fig. 7. 68% and 95% c.l. constraints for the joint variation 
of a and m e from different data sets. 
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Fig. 6. 68% and 95% c.l. constraints for m e /m e o versus 
fifc/i 2 , for the analysis of the variation of m e alone. Results 
from different data sets. 




mass m e during recombination epoch, using the latest CMB 
data, including the 7-yr release of WMAP. We perform sev- 
eral statistical anal yses adding two d ifferent data sets; the 
H prior taken from lRiess et all (|2009h : and the power spec- 
trum from Sloan Digital Sky Survery DR7 LRG (|Reid et all 
I2009D . The bounds on the variation of the constants are 
tighter than previous results because of the higher preci- 
sion of the new data used in this work. 

Our results show no variation of the constants at re- 
combination time. We also emphasize that the constraints 
depend strongly on which data set we choose in the analy- 
sis, due to the large degeneracy between a or m e and Hq. 
Yet, the results are consistent within 1-a. 
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mean value of m e is more affected than the mean value of 
a. These results can also be seen in Table [5J 

A vari ation of a or m e affec ts the recombination sce- 
nario (see IScoccola et al.l (|2009[ ) for example). As a con- 
sequence, the angular diameter distance at recombination 
is modified if any of these constants varies. This results 
in a change in the Dopp ler peak positions and heights (see 
Kaplinghat et al. (1999) for example). This explains the de- 
generacy between a and m e shown in Fig. [7] and confirmed 
by the correlation coefficient. On the other hand, the de- 
generacy between a or m e with the baryon mass density 
or the Hubble constant can be explained since these effects 
are similar to a change in the cosmological parameters. 



4. Conclusions 



In this paper we have updated the constraints on the time 
variation of the fine structure constant a and the electron 
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